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Abstract 

We study partial rate asymmetry in decays with single isospin final state. For 



K meson or hyperon decays, the partial rate asymmetries are always zero if the 
final states are single isospin states. In B decays the situation is dramatically 
different and partial rate asymmetries can be non-zero even if the final states 
are single isospin states. We calculated partial rate asymmetries for several 
B decays with single isospin amplitude in the final states using factorization 
approximation. We find that in some cases the asymmetries can be large. 
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Typeset using REVTpjX 



CP violation is one of the few remaining unresolved mysteries in particle physics. The 
explanation in the Standard Model (SM) based on Cabibbo-Kobayashi-Maskawa (CKM) 
matrix JIJ is still not established, although there is no conflict between the observation of 
CP violation in the neutral K-system and theory |2[. It is important to carry out more 
experiments to test CP violation in the SM. The study of CP violation in the B system 
is very important which may provide crucial information about CP violation ||. The B 
system offers several final states that provide a rich source for the study of this phenomena. 
In many cases CP violation in B decays occurs in a quite different form from K or hyperon 
decays. In this paper we study CP violating partial rate asymmetries in B decays. We clarify 
some subtleties for CP violation in partial rate asymmetry in relation to isospin analysis. 

The decay amplitudes for Kaon or hyperon are customarilly parametrized according to 
isospin decomposition in the final states because isospin states are eigenstates of strong 
interaction. For example the amplitude for K° — > 7r + 7r~ decay can be parametrized as 

A = A e iS °- +5 ° + A 2 e l5l+S " , (1) 

where superscript and 2 indicate the isospin of the final states. 5 l w and 5 l s are the CP 

violating weak and CP conserving strong phases, respectively. 

The partial rate asymmetry A asy is given by 

|A| 2 -|A| 2 = 2AAsin(^ ~ <%)sin(ff - £ 2 ) 

asy \A\ 2 + |A| 2 " Al + A\ + 2A^ 2 sin(($o - S 2 W )sin(5° - 5 2 s ) ' 1 ) 

The conditions for non-zero partial rate asymmetry in this case are: there must exist at least 

two different isospin decay amplitudes with different weak, and strong rescattering phases. 

It is clear that for final states with single isospin, for example K~ — > 7r~7r°, the partial 

rate asymmetry vanishes. The same equation could be used for B — > txtx decays. But the 

argument is now incorrect. We wish to consider the difference in some detail. 

We now consider the same process as above from the quark level. In the SM the effective 

Hamiltonian responsible for K — > 7r + 7r~ at the quark level can be parametrized as ||, 

Heff = £[K S + v cs v; d ct + VtsV^o, , (3) 



where q are Wilson Coefficients (WC) of the corresponding quark operators Oj. At the 
one loop level c" contain the tree and u internal quark contributions, q'* contain internal 
c- and t- quark contributions. Since the u, u pair is lighter than s-quark, at the one loop 
level with u quark in the loop, absorptive amplitudes will be generated, and cf has a strong 
rescattering phase |J. On the other hand no absorptive parts exist for cf f . Naively, one 
would obtain partial rate asymmetry for K° — > 7r + 7r _ which seems to have nothing to do 
with strong phase shifts in different isospin amplitudes. However, this argument turns out 
to be wrong. It has been pointed out by Gerard and Hou || that CPT theorem is violated 
if one is not careful to include all diagrams of the same order 0. The interference term 
responsible for rate asymmetry due to c" is an interference between penguin amplitudes of 
order a 2 s . To this order there is also another contribution which is the interference between 
the tree amplitude and the higher order penguin amplitudes with absorptive part developed 
in the vacuum polarization of the virtual gluon. This contribution cancels the previous 
contribution. In practical calculation, cf must be treated as real in this case. The rule is 
that for the phase of any one of the penguin WC, if there is a tree amplitude with the same 
CKM factor, the phase in that penguin WC must be removed when the final states are 
the same as the tree amplitude. In a more general formulation of the problem from CPT 
theorem and unitarity considerations, Wolfenstein showed that any diagonal strong phase 
(the phase due to rescattering of the states which are the same as the final states) do not 
contribute to partial rate asymmetry ||. The phases in cf are diaganal phases in the present 
case ||. 

At the hadron level, the decay amplitude is given by 

A{K v _> 7r+7r - ) = J2< n + n-\[V: d V us c« + V^V^ + V£y u <$\Oi\K > . (4) 

i 

The strong rescattering phases are generated by rescattering the two pions in the final states. 
This is because in this case, only two pion final states are opened in this kinematic region 
with the right parity. Three pions in the intermediate states are allowed kinematically, but 
three pions would not rescatter into two pions because of G parity conservation. The part 

3 



responsible for absorptive amplitude is given by 

A{K o _> ff+7r -) = J2J2< ^ItO^ + V; d V cs c1 + VZVutftlimr)! > 

i I 

x < (wMOilK > , (5) 

where I is summed over isospin eigenstates. 

Since isospin symmetry is respected by strong interaction, we can generally parametrize 
the hadronic matrix elements as, after the rescattering strong phase shifts are included, 

£j < * + K-\[V^V us c» + V; d V cs ct + V t * d Vt s cl]\(iTir)i >< faUO^K > 
= KM + V: d V cs 4 + V td V ts cj} [x%e i5 ° + x 2 a 2 e i5 °] , (6) 

where a 0t2 and 5®' 2 are the isospin eigen-amplitudes and strong rescattering phases from all 
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contributions, and x{ are the Clebsch-Gordan coefficients for each operator. Knowing that 
the absorptive parts in cf are canceled by the other effects, we should take all e"' c '* in the 
above equation to be real. We have for the isospin amplitudes A 02 

A) = EtK^Kscr + v; d v cs ct + V td V ts 4]x%e l5 ° , 

i 

M = EtO^ + V: d V cs c\ + V td V t 4]x 2 t a 2 e^ . (7) 

i 

Since in general x® is not equal or proportional to xf, Aq and A 2 do not have the same weak 
phases 

SS = Arg{Y\V: d V us dt + V^V^ + V td V ts c^) , 

i 

6™ = Arg(J2[V: d V us ^ + V^c* + V^V^x 2 ) . (8) 

i 

Thus the decay amplitude for K° — > -k + -k~ can be parametrized in the form in eq. ([!]). The 
discussion can be easily generalized to hyperon decays. 

For K~ — > 7T _ 7r°, the final state has only 1 = 2 amplitude, the decay amplitude is of the 
form, 

M = EtKdKscr + V; d V cs ct + V td V ts 4]x 2 ~a 2 e l5 ° . (9) 

i 
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It is clear that the particle-antiparticle rate asymmetry vanishes. 

For B — > 7T7T the situation is, however, very different. We will use the same notation 
for the effective Hamiltonian. Of course we should keep in mind that now the operators 0« 
contain a b quark. Now because the b quark is heavier than a u, u pair, and also a c, c pair, 
both c", c\ have strong rescattering phases at the one loop level. At the hadron level, the 
part responsible for absorptive amplitude is given by 

a(b° - Tr+Tr-) = EE < « + *-\[v: d v ub $ + v: d v cb ct + v; d v tb 4]\i >< iio^b > 

i I 

+ EE < TT + TT-IO^I/c >< Ic\O t \B° > , (10) 

i I c 

where I is summed over non-charmed on-shell intermediate states like, n + 7i~ , 7r°7r°, and 
I c is summed over charmed on-shell particle intermediate states like, DD etc. In this case 
we can remove the phases in c" because in this decay there is a tree amplitude with the 
same CKM factor. However, now there is another class of phase shift due to < i\ Jr i\~\c c i \I c > 
which can not be removed. The phase shift of this type to the lowest order is due to the 
absorptive part in c\. Because these new phases are generated by charmed particles in the 
intermediate state, these phases will only appear in the term proportional to V cb V* d . We 
obtain the / = 0, 2 decay amplitudes A' 2 for B — > nn 

A'o = YSVudVub^ + V: d V ch c\ + V t * d V tb c*]x? a e^° , 

i 

A' 2 = EKdK fe cr + V: d V cb ct + V td V tb cl\x? a 2 e< 2 . (11) 

i 

Now we should treat c"'* to be real and only c\ to be complex (non-zero rescattering phases) 
up to order a 2 s in the asymmetry. 

Let us now compare these amplitudes with the amplitudes for K — > mr. First we note 
that because in the case for B — > nn more on-shell intermediate states are allowed, i.e. nn, 
-kit-kit etc, the rescattering phases 5' s ' 2 include inelastic channels unlike the elastic phase 
shifts 5®' 2 . Second we note that for B — > ixix there are additional strong rescattering phases 
due to on-shell charmed intermediate states. The strong rescattering phases in c\ are not 
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canceled by any other contributions. This has a very important consequence that particle- 
antiparticle rate asymmetry can occur in a single isospin amplitude. In fact this happens 
quite often in B decays ||. In the following we study four different processes representing 
different types of B decays, b — > tps, b — > <ps, B~ — ► ^7r _ , and B~ — ► i]K~ . 

When considering partial rate asymmetry for single isospin final state in B decays, we 
do not need to know the values for the overall strong rescattering phases. We, however, 
need to know the relative strong rescattering phase shifts between amplitudes with different 
CKM factors by calculating various on-shell rescattering processes. This calculation is very 
difficult to carry out. However, we believe that the WC's and the phases calculated at the 
quark level could be good indications of the sizes and the signs of the strong rescattering 
phases by appealing to duality. The absorptive parts of hadronic processes are given quite 
accurately by considering the corresponding quark loops as in the calculation of R in e + e~ 
scattering. In our later calculation, we will use this approximation. 

In the SM the amplitudes for B decays are generated by the following effective Hamilto- 
nian: 

Q 10 

Hlff = ^[VfbVf^ctOy + c 2 O q 2f ) - YSYubV^ + V ch V: q c\ + V tb V*cj)Ol} + H.C. , (12) 

where the superscripts u, c, t indicate the internal quarks, / can be u or c quark, q can be 
d or s quark depending on if the decay is a AS = or AS = — 1 process. The operators Of 
are defined as 

O q fl = q a j ti LupU /3 j> J 'Lb a , Of/ = qj^Luu^Lb , 

01,5 = Ql^bq'^L(R)q' , Of j6 = q a -f^Lb p q'^^L(R)q' a , (13) 

3 3 
Ot,9 = ^Q7^ L be q >q'j'"R{L)q' , O| 10 = -q^^Lbpe^q'^^R^q^ , 

where R(L) = 1 + (— )7s, and q 1 is summed over u, d, and s. 0\ are the tree level and QCD 
corrected operators. O3-Q are the strong gluon induced penguin operators, and operators 
O7_io are due to 7 and Z exchange, and "box" diagrams at loop level. The WC's c{ are 
defined at the scale of /1 ~ m& which have been evaluated to the next-to-leading order in 
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QCD Hp, 10]. We give the non-zero coefficients below for m t = 176 GeV, a s (mz) = 0.117, 



and n = m b = 5 GeV, 

d = -0.307 , c 2 = 1.147 , 4 = 0.017 , c\ = -0.037 , c\ = 0.010 , c\ = -0.045 , 
4 = -1.24 x 10~ 5 , 4 = 3 - 77 x 10 ~ 4 > 4 = -0.010 , 4o = 2- 06 x 10 ~ 3 » 
$$ = -<Z$/N = P:/N, c u 7 $ = Pr, cZo = (14) 

where A 7 is the number of color, c* are the regularization scheme independent WC's ob- 



tained in Ref. fllP] . The leading contributions to P* e are given by: Pj = (a s /87r)c 2 (10/9 + 
G(mi, n,q 2 )) and P* = (a em /97r)(A^ci + c 2 )(10/9 + G(rrii, /x, g 2 )). The function G(m,fi,q 2 ) 
is give by 

^/ f 1 , . m 2 — x(l — x)q 2 , , 

G(m,fi,q 2 ) =4 / x(l -x)dxln ^- . (15) 

jo iJtr 

All the above coefficients are obtained up to one loop order in electroweak interactions. 
When q 2 > 4m 2 , G(m, fi,q 2 ) becomes imaginary. In our calculation, we will use m u = 5 
MeV, rrid = 10 MeV, m s = 175 MeV, m c = 1.35 GeV, and the averaged value ml/2 for q 2 . 

We must be careful in using absorptive parts of the above WC's. We should always 
remove phase shift discussed previously, otherwise one would obtain results violating CPT 
theorem. 

To obtain exclusive decay amplitudes, we need to calculate relevant hadronic matrix 
elements. Since no reliable calculational tool exists for two body modes, we shall use factor- 
ization approximation to get an idea of the size of asymmetry A asy . The numerical numbers 
obtained should be viewed as an order of magnitude estimates. The important message is 
that CP violating partial rate asymmetry in some single isospin channel decays are indeed 
non-zero and can reach significant magnitude. 
Partial rate asymmetry in b — > ips 

Since ip carries no isospin, the final state is a single isospin state in this decay. At the 
hadronic level, this decay includes B — > ipK, B ipK*, and etc. This decay is particularly 
interesting because it has a large branching ratio. CP violating partial rate asymmetry in 



this type of decay was first studied in the early 80's by Brown, Pakvasa and Tuan ||11|| . Let 

us now analyze this asymmetry using factorization approximation. In this approximation, 

the partial rate asymmetries for b — > ips and B~ — > ipK~ are the same. We have 

G c 
A{b^ ^s) = ^=*7m(1 - 7s)H< ^\u^u\0 > V ub V: s ( Cl + ^) 

+ < i>\*fc\o > [v cb v; s ( Cl + 1) - E + § + 4 + 1 + 4 + § + 4 + f )]} ■ (ie) 

Here the first term corresponds to an annihilation contribution which is usually small. If 
this term is neglected, we have 

G 

A{b -> if;s) = 7^«7^(1 - 7s)& < V|57 M c|0 > 

x [VkC(C7 - - c«) - V cb V; s (C + d + | - c? - eg)] ; (17) 

where C = c| + c| + + c\ + (c| + Cg + c\ + c\ )/N. Any phase shift in should be removed 
from our previous discussions. Only the absorptive part in c" 9 generate effective strong 
rescattering phases. The interference which cause the partial rate difference in this case is 
of order a em instead naively expected a s because the strong penguin generated absorptive 
amplitude cancel in Cg 5 + c^q/N. In our numerical calculations we will use N = 2 favored 
by experimental data, and \V US \ = 0.2205, |V c b| = 0.04, and \V u b/V c b\ = 0.08. We find that 
the partial rate asymmetry is less than 10 -4 . The result is shown in Figure 1. One can 
easily obtain the asymmetry for b — ► tpd by scaling the asymmetry by a factor of \V c b/V u b\ 2 . 
The asymmetry in this case is much large but the branching ratio is much smaller. In Ref. 
12fl using absorptive amplitudes generated by rescattering color octet state, asymmetry was 



estimated for b — > ipd. Using the same calculation for b — > ips, we find that the partial 
rate asymmetry is about the same order of magnitude obtained here. This asymmetry is 



an order of magnitude smaller than that obtained in Ref. OJT3I]. Note that here we have 



neglected the annihilation contribution estimated in Ref. flTTffSf. However, it is found that 



when current knowledge about CKM parameters is used, the caculation of Ref. [ |TT1 , |I"3f also 
yields an asymmetry below 10~ 4 . Hence including the annihilation diagram is not going to 
change our result significantly. 
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Partial rate asymmetry in b — > <ps 

The final state is a single isospin state because <fr is isospin singlet. This decay is induced 



by pure penguin interaction. CP violation in this process was first evaluated in Ref. |T3 
This process is not affected by the previously mentioned effect. The partial rate asymmetry 
is generated by interference of different penguin amplitude which is of order a 2 s . We have 

G 

A(b -> 0s) = ^=^V(1 - 75)6 < 0|s 7 ^s|O > 

x [V u bV* s (C - 2(c« + |) + 4(1 + 1) + <%) 

+ V cb V; s (C - 2(4 + I + 4(1 + 1) + eg)] , (18) 

where C" = 2(4 + c| + c|) + 2(4 + c| + c\)/N - (4 + 4 + 4o + (4 + 4 + c \o)l N )- In this 
case we find the partial rate asymmetry is of order O(10~ 3 ). The result is shown in Figure 
2. 

Partial rate asymmetry in B~ — > ?]7r 

This is an exclusive decay. Because 77 is an J = particle, the final state is a single 
isospin state with 1 = 1. In the factorization approximation, we obtain 

A(B- -> ,770 = ^[V^Vi^CT- + T^ v D ut ) 

+ V cb V: d (Tl n C ct + Tl v D ct )}, (19) 



where 



^- 1 + |-f-4- + |(4 + f) + j(| + 4o) 

+ ^(-2|-24 + 4 + |-4-^) + {4--4}, 
^ = | + c 2 -(| + 4 + | + 4o + 2^(| + 4 + | + 4)) + {4--4}, (20) 

where = mi/ (2md(mb — m^)) and X v = m\j{m u + m,j){mb — ma). C ct and D ct are 
obtained by setting ci j2 to be zero and replacing the superscript u by c. T^. and Tg are 
defined as, 

TjBir = < TT-ldj^l - l5 )b\B- >< ~ 75)«|0 > 



= l (fr,/V3)Ft(m 2 v )(m 2 B -ml), 
Tb v = < - n)b\B- < n-d^(l - 75 )u|0 > 

= i(fjV3)F B «(ml)(m 2 B -m 2 ri ), (21) 

where F Bp (q 2 ) are the transition form factors between B meson andp meson (where p could 



be 7r or 77) defined in Ref. [p~5|,il6|] , f n = 93MeV, and we will use f w ~ f v . This time any strong 
phase in c" must be removed. The interference causing partial rate difference is of order a s . 
The result is shown in Figure 3. In the numerical calculation, we have included the 77 — rf 
mixing effect with the mixing angle 9 = —20°. The figures are ploted with f m = f m = f w . 
The results are not sensitive the the mixing effect. The asymmetry can be quite large. The 
branching ratio is 3 to 4 times larger (O(10~ 5 )) if the form factors in Ref. [15] are used than 
the one obtained using the form factors in Ref. [[K . 



Partial rate asymmetry in B — ► rjK 

The final state is a pure 1 = 1/2 state. We have 

A(B- -> rjK-) = ^[V ub V: s (n K C ut + 



u u o u u 

N N 4 N wJ 2 K 7 N 9 N ' 

„u u 

+ 2X.{-± + 4 - -± - cj) + {c* -cj} 



B "' = l + C2 -|- c:; - 2 ^ ( l + < + F + C " 

c, 



- ^ - c? + R - -c<} , (22) 

where = m 2 K /(m s + m u ){m h - m u ), X s = m 2 /(2m s (m b - m 8 )), T% K = 

i(f v /V3)F BK (m 2 )(m 2 B - m 2 K ), = i(f K /^)F Br, (m 2 K )(ml - m 2 ). Similarly the coef- 
ficients C ct and D ct are obtained by setting Ci )2 to be zero, and replacing the superscript u 
by c in C ut and D ut . This time the the interference term causing partial rate difference is, 
again, of order a s . The result is shown in Figure 4. In this case the results are sensitive to 
the mixing effect. Within the allowed ranges for fn^^g, the branching ratio can change by 
a factor of 5. The branching ratio can be as large as 4 x 10 -6 using the form factors in Ref. 
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15 1, and it is smaller by a factor of 3 to 4 using the form factors in Ref. [[Uj. The mixing 



effect on the asymmetry is less sensitive. Again, the asymmetry can be quite large. For the 
same values of form factors, we agree with the results obtained by Du and Guo in Ref. || . If 
the electroweak penguin effect is neglected, we also agree with Kramer, Palmer and Simma 
in Ref. |J when the same form factors are used. 

To conclude, we have shown that partial rate asymmetry in decays of particles containing 
a strange quark with single isospin final states, is always zero because the allowed intermedi- 
ate on-shell states are limited. However, the situation in B decays is dramatically different. 
In the latter case, more intermediate on-shell states with different CKM factors are allowed, 
CP violating partial rate asymmetries need not to be zero even if the final state contains 
only a single isospin state. 

We have carried out detailed analyses for four types of B decays. The partial rate 
asymmetry in b — > ips is small because the interference causing rate deference in particle 
and anti-particle decay rates is of order a em due to cancellation. The partial rate asymmetry 
in b — > <ps is also small (O(10 -3 ). The asymmetry in exclusive decays considered here are 
larger. The partial rate asymmetry in B~ — > rj-K^ and B~ — > r]K~ can be quite large. 

This work was supported in part by the Department of Energy Grant No. DE-FG06- 
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Laboratory for hospitality where part of this work was performed. 



11 



REFERENCES 

[1] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. 
Theor. Phys. 49, 652 (1973). 

[2] A. Buras, Preprint, MPI-Pht/95-30, TUM-T3 1-88/95, hep-9504269; J. Rosner, 
Preprint, EFI-95-36, hep-9506364; A. Ali and D. London, DESY 95-148, UDEM-GPP- 
TH-95-32, hep-9508272; 

[3] For a review see: B Decays, edited by S. Stone, World Scientific, 1994. 

[4] M. Lautenbacher and P. Weisz, Nucl. Phys. B400, 37(1993); A. Buras, M. Jamin and 
M. Lautenbacher, ibid, 75(1993); M. Ciuchini, E. Franco, G. Martinelli and L. Reina, 
Nucl. Phys. B415, 403(1994) 

[5] M. Bander, D. Silverman and A. Soni, Phys. Rev. Lett. 43, 242(1979). 

[6] J. Gerard and W.-S. Hou, Phys. Rev. Lett. 62, 855(1989). 

[7] J. Gerard and W.-S. Hou, Phys. Rev. D43, 2909(1991); H. Simma, G. Eilam and D. 
Wyler, Nucl. Phys. B352, 367(1991). 

[8] L. Wolfenstein, Phys. Rev. D43, 151(1990). 

[9] G. Kramer, W. Palmer and H. Simma, Z. Phys. C66, 429(1995); D.-S. Du and L. Guo, 
Preprint, BIHEP-TH-95-28. 

[10] N. G. Deshpande and Xiao-Gang He, Phys. Lett. B336, (1994)471. 

[11] T. Brown, S. Pakvasa and S. F. Tuan, Phys. Lett. 136B, 117(1984). 

[12] J. Soares, Phys. Rev. D52, 242(1995). 

[13] L.-L. Chau and H.-Y. Cheng, Phys. Rev. Lett. 53, 1037(1984). 
[14] N.G. Deshpande and J. Trampetic, Phys. Rev. D41, 2926(1990). 

[15] A. Deandrea, N. di Bartolomeo, R. Gatto and G. Nardulli, Phys. Lett. B318, 549(1993). 

12 



M. Bauer, B. Stech and M. Wirbel, Z. Phys. C34, 103(1987). 



13 



FIGURES 
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FIG. 1. The partial rate asymmetry for b — > tps. The vertical axis is the asymmetry and the 
horizontal axis is the value in degree for the phase angle 7 in the Wolfenstein parametrization. 
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FIG. 2. The partial rate asymmetry for b — ► 4>s. 
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FIG. 3. The partial rate asymmetry for B — > r/ir . 
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FIG. 4. The partial rate asymmetry for B — > r/K . 
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